Spin resolution of the electron-gas correlation energy: 
Positive same-spin contribution 
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The negative correlation energy e c (r s , £) per particle of a uniform electron gas of density parameter 
r 3 and spin polarization £ is well known, but its spin resolution into H, ||, and J. J. contributions is 
not. Widely-used estimates are incorrect, and hamper the development of reliable density functionals 
and pair distribution functions. For the spin resolution, we present interpolations between high- and 
low-density limits that agree with available Quantum Monte Carlo data. In the low-density limit 
for £ = 0, we find that the same-spin correlation energy is unexpectedly positive, and we explain 
why. We also estimate the f and J. contributions to the kinetic energy of correlation. 



The uniform electron gas is a paradigm for density 
functional theoryji* 2 ^ the most widely-used method for 
electronic structure calculations in both condensed mat- 
ter physics and quantum chemistry. The effects of ex- 
change and correlation can be evaluated and under- 
stood in the uniform-density limit, and then transferred 
to realistic systems. This is done not only in the lo- 
cal spin density (LSD) approximation but also beyond 
LSD in generalized gradient approximations (GGA's), 
meta-GGA's, and hybrid functionals £ The correlation 
energy e c (r s ,£) per particle in a uniform gas of den- 
sity parameter r s — (47maQ/3) -1 / 3 and spin polarization 
£ = (n-f — n\,)/n (where n a is the density of spin-er elec- 
trons and n = n-f + ) is well known, for example from 
Quantum Monte Carlo (QMC) studies 4,5 that have been 
parametrized£iL& to respect known limits, but the spin 
resolution of e c into tj, ft, and J, J, contributions is not 
known. In this work, we determine the spin resolution 
for all r s and £ as an interpolation between high- and 
low-density limits, consistent with ( = QMC data£ 

This spin resolution is of interest in its own right, 
and can also be used in several ways: (i) Some beyond- 
LSD correlation energy functionals need a missing spin 
resolution^ or have been constructecUSiiiiiSii^ on the ba- 
sis of the exchange-like ansatz of Stoll et a/jlA 

[n h nj] « E c [n T , - E D [n h 0] - E c [0, , ( 1 ) 

for the uniform gas. This assumption was shown (using 
QMC results) to be inaccurate for ( — (see Fig. 
in Ref. flH although the significance of this observation 
for density functional theory was not fully recognized 
there. Our work provides a firmer basis than Eq. 
for such constructions, (ii) Correlation energy function- 
als such as the local spin densitji and generalized gra- 
dient approximations^ etcH can alternatively be con- 
structed without a spin resolution, but their later spin 
resolution (to permit comparison or combination with 
correlated-wavefunction resultsiS* 1 ^* 1 ^) demands such a 
resolution for uniform densities, (iii) A sophisticated an- 
alytic model^I is now available for the pair distribution 
functionl^fliSi g xc (r s , £, it) of the uniform gas for all r s 



and £. Our present work provides the missing ingredi- 
ent needed to find the corresponding spin-resolved pair 
distribution function, which could serve as the starting 
point for the development of density functionals such as 
spin-resolved weighted density approximations.— (iv) An 
estimate can be made for the £ dependence of the f and 
I contributions to the kinetic energy of correlation, a key 
ingredient for the approach to spin dynamics of Qian and 
Vignale— and also for the momentum distribution 2 ^ of a 
spin-polarized electron gas. 

We shall first derive exact limits at high densities 
(r s — > 0) and extreme low densities (r s — > oo) using sim- 
ple physical arguments. In the latter limit, we find that 
the same-spin contribution to the correlation energy can 
be positive, and we provide an intuitive physical picture 
to explain this feature. While the total correlation en- 
ergy must be negative, individual terms of it (e.g., the 
kinetic energy of correlation) can be positive. We then 
build up and discuss our interpolation formulas. 

Definitions - Correlation effects arise from the 
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FIG. 1: Fraction of || correlation energy, Ffi(r s ,Q = 
eP(r s , (")/e c (r s , f) at £ = 0. Our Eq. © is compared with 
the GSB-^ values extracted from QMO^ data (•), and with 



the Stoll et al^ PW9&&, and SKTW 
lence electrons have 2 < r s < 6. 
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FIG. 2: Spin resolution F aa i(r 3 , f) = e" CT (r s , C)/e c (r s , C) as a function of £ for different r a . The high-density (HD) and low- 
density (LD) limits are given in Eqs. © and @. The r 3 = 3.28 curves correspond to the SKTB-& scaling relation of Eq. |0J , 
while for other density values (r s — 1, 10 and 100) our interpolation formulas of Eq. @ have been used. 



Coulomb interaction, which is a two-body operator. 
When evaluating the energy of the system, (\f'|i/|\f'), one 
can split the sum over the electron spins into ||, II, an d 
|| contributions. The corresponding splitting of the cor- 
relation energy of the uniform electron gas, 



Jit 



(r s X) + 4HrsX) + 4 l (rsX), (2) 



1 1 / 



is the object of this paper. The real-space analysis of 
the spin-resolved correlation energies e° a (r s ,C) is pro- 
vided by the correlation holes n a >'g^' T (r s ,£,u) (see, e.g., 
Ref . l2l|) . where u — |ri — r2 1 is the electron-electron dis- 
tance: 



(r s ,C) = 27r^ 
n 



n<r']% C (r s X,u)udu. (3) 



The correlation hole n a i~g°' y (r s , £, u) describes the 
change (due to correlation only) of spin-tr' electron den- 
sity at u, when a spin-cr electron is at the origin. g c is 
averaged over coupling strength, while g c is for full cou- 
pling strength. We define fractions F aa i (r s , () such that 



(r s ,() = e c (r s ,0F r7r7 ,(r s ,0, 



(4) 



and we investigate their properties. In what follows, 
we use Hartree atomic units, and the parametrization 
of e c (r s , C) and its limits from Ref. 0. 

Exact limits - When r s — > 0, the Coulomb electron- 
electron interaction can be treated as a perturbation to 
the non-interacting Fermi gas. The first-order (in the 
Coulomb potential) correction term gives the exchange 



energy e x = 
Al = 3_ 



+ e^, where 4T = -^(1 + C) 4 / 3 , 

, , ,_„;(! -C) 4/3 , and a = (97T/4)" 1 / 3 . As for cor- 
relation, the real-space analysis of the exchange energies 
is provided by the exchange holes n rr \qZ a (C : u/r„) — 1], 
which are analytically known (see, e.g., Ref . 12111 . 

The second order correction to the energy of the non- 
interacting Fermi gas is the sum of a direct term and a 
second-order exchange term. Only the direct term di- 
verges, and, when a cutoff cx (due to Thomas- 
Fermi screening effects) at small wavevectors is intro- 
duced, gives rise to a leading term in e c (r s , () , equal to 



co(C)hir s . The function cq(C) is exactly known^ The 
direct term (Eq. (5.110) of Ref. l26|) can be divided into 
||, ||, and || excitation pairs to derive 
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pHD 
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i + C 

47(C)' 



(•5) 



with 7(0 = Co(0/cq(0), as conjectured in Ref. 25. (Since 
Fu(r s ,0 = F n( r s,~C) andF n = 1-F n -F u , we only 
report formulas for ||.) The Stoll et al. ansatz of Eq. JJJ 
is thus correct for r s — > (and for all r s when |£| = 1, 
but not otherwise). 

In the opposite or strong-interaction limit, r s — > oo, 
the long-range Coulomb repulsion between the electrons 
becomes dominant with respect to the kinetic energy, 
and thus with respect to statistics; Coulomb repulsion 
suppresses electron-electron overlap so that the electrons 
no longer know they are fermions. In this limit, the to- 
tal energy becomes independent 4,5,8,21 - 27 of £. Its lead- 
ing term in the r s — > oo expansion is equal to —d\/r s , 
where& d\ ~ 0.892, and is purely potential energy, with 
no kinetic energy contribution. In this limit, the total 
energy is thus equal to the exchange-correlation energy 
£xc = £x + £c- Moreover, since the statistics becomes 
irrelevant, we expect that 



i + C 



r [ [ = 
xc 



i-C 



(i-C 2 ) x 



f Tl = e 



(6) 



where the prefactors take into account the available 
numbers of pairs. In other words, we expect that 
/ °° du 4:iru 2 g°c l u becomes independent of a and a', 
so that spin structure becomes unimportant for the 
exchange-correlation and total energies (although very 
important for the correlation energy alone). Then the 
Foo> (r s -» oo, = Faa'iO ar e given by 



F LB (n = 3(1 + O 4 / 3 - 27ra(l + 

TT lU 3[(i + 4/3 + (i-0 4/3 ] 



2 di 



8na d\ 



(7) 



The high- and low-density F aa i are displayed in Fig. [21 
We see that, in the spin-unpolarized gas, the same-spin 
(11+11) contribution to the correlation energy is 50% 
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FIG. 3: Upper panel: the spin-resolved pair distribution func- 
tions for the paramagnetic gas. The dashed arrows show the 
trend of the holes as the coupling strength r a is increased. 
Lower panel: real-space analysis of the correlation energy in 
the extreme low-density limit, for the paramagnetic gas. The 
results are from the model of Ref. l2ll 



when r s — > but roughly when r s — > oo. This can 
be understood in a simple way. The exchange hole seen 
by the same-spin electrons is deep for electron-electron 
distances u < r Sl as shown in the upper panel of Fig. [3] 
(solid line, + J, J,). But there is a second length scale, 
the Thomas-Fermi screening length ^/r^. For r s — > 0, 
the important correlations, which determine the leading 
term (cx lnr s ) of e c , arise from this second length scale, 
y/FH 3> r s , and are essentially unaffected by exchange: 
the electrons that participate in this correlation have no 
way to know if the electron at u = is spin-| or spin- J,, so 
by symmetry the same-spin and opposite-spin correlation 
energies are equal. In the opposite limit, r s — ► oo, the 
antiparallcl-spin correlation hole can get deep for u < r s , 
as shown in the upper panel of Fig. [3J 

As r s increases, g\\. deviates more and more from its 
non-interacting value (equal to 1 for all u), the only con- 
straint being its positivencss. But the same-spin corre- 
lation hole is "blocked" from doing this by the exchange 
hole (see, again, the upper panel of Fig. |3J). Thus the 
system minimizes its energy by focussing the correlation 
on opposite-spin pairs. In the extreme low-density limit, 
a simple qualitative picture can be obtained by using the 
correlation-hole model of Ref. |^] (in which energetically 
unimportant long-range oscillations are averaged out); 
in the lower panel of Fig- EI we report the corresponding 
real-space analysis of e]} and eJ T + for r s — > oo. We 
see that the same-spin correlation hole for u <r s cannot 
get as deep as the opposite-spin one. 

Figure|21also shows that in the spin-unpolarized gas the 
same-spin correlation energy is slightly positive (F aa < 
0) when r s — > oo. In this limit, the electrons correlate 



strongly, and the exchange-correlation holes show a high 
first-neighbor peak at u « 2r s (lower panel of Fig. |3J). 
If the only effect of same-spin correlation were to push 
same-spin electrons away from the region of small u and 
pile them up at u « 2r s , then (by the sum rule integral 
J duAiru 2 n a g°' J (u) = 0) the same-spin correlation energy 
[Eq. © with <r — a'] would necessarily be negative. So, 
what must really happen is that the same-spin electrons 
that accumulate in the peak at it « 2r s include some that 
have been pushed out from u <§; 2r s and some that have 
been pulled in from u 3> 2r s . This is again illustrated 
in the lower panel of Fig. |3J We interpret the second 
zero of g c , which appears at large u but only at large r s , 
as the energetically important remnant of the long-range 
oscillation of g c in a Wigner crystal. 

Positive same-spin correlation energy may be an ex- 
otic effect, but the blockage of negative same-spin corre- 
lation also occurs in a non-magnetic Mott insulator, e.g., 
an expanded lattice of hydrogen atoms where Coulomb 
correlation suppresses the (Is) 2 configuration on a given 
site. The blockage of same-spin correlation occurs even 
in a weakly-correlated spin-unpolarized system when 
the correlation hole is spatially constrained, as for an 
a ^ om> i2 i i8 i i9 j n neon a tom, the true (as cited in 
Ref. Il8|) anti-parallel-spin correlation energy is 65% of 
its LSD value, while the true parallel-spin correlation en- 
ergy is only 30% of its LSD value. 

Interpolation between high and low density - We want 
to build up interpolation formulas for F a<y i (r s , £) that 
include all the information available on the spin resolu- 
tion of e c . Besides the high- and low-density limits, we 
have data for F arT i(r s ,0), in the range 0.8 < r s < 10. 
These data have been obtained in Ref. ^| (GSB) by in- 
tegrating spin-resolved QMC correlation holes£ More- 
over, Schmidt et alw& (SKTP), starting from nearly-exact 
limits of the spin-resolved correlation holes, proposed a 
scaling relation that is in agreement with the GSB data 
at r s = 3.28, and that, as shown in Fig. [21 (curves la- 
belled with "r s = 3.28"), lies in between the high- and 
the low-density limits with a very "reasonable" shape. 
The SKTP scaling should thus be a good "intermediate 
point" for our interpolation formulas. We thus define 
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(0 = 



i + C 



11/6 



,(3-28,1) 
=(3.28, 0' 



(8) 



and we parametrize F aa i (r s , £) as 



1 + C y/r s + B r s 



Aaa>(() is found by requiring that F aa < (3.28, C) 
^ TP (C), i.e., 



(9) 



vW(C) 



— — +c (c) + 

B VX28 [^P(() - F%(Q]. (10) 

The form of Eq. @ is motivated by the expression for the 
correlation energy given in Ref.|7|. The parameters B and 



4 



C are fixed by a best fit of F a<J i (r s , 0) to the GSB data for 
r s £ [0.8,10]: B = 0.178488, C = 2.856. In Fig. □ our 
F||(r s ,0) is compared with the GSB data^ and with 
the widely- used Stoll et alM* ansatz of Eq. (JIJ, which 
strongly underestimates the fraction of || correlation en- 
ergy at metallic and lower densities. The results for the 
paramagnetic gas corresponding to other proposed scal- 
ing relations are also shown. Our interpolation formulas 
as functions of C, at r s — 1, 10, and 100, are displayed 
in Fig. H 

Kinetic energy of correlation - Defining 29 ej = ejj + 
ie^ (with a similar equation for |), the adiabatic con- 
nection between the non-interacting and interacting lim- 
its for a given density suggests estimating the f and j 
contributions (from the one-particle density matrix) to 
the kinetic energy of correlation t c = t\ + t\. asS 

^(r.,C)»-^M(r„C)], (11) 

although as Ref . |3(3 points out there is only one coupling 
constant with a Hellmann-Feynman theorem, not one for 
each a. Taking Eq. (|ll|l as a plausible approximation, we 
find that the corresponding result for t\ — 1\ is in reason- 
able agreement with the scaling relation given in Eq. (29) 



of Ref. [2Ej. (For r s < 5, the difference is less than 3.5%). 
Via Eq. (|ll|l . we also confirm that, for 1 < r s < 10, the 
quantity (t\ —t*)/t c (r s , C) is almost independent of r s , as 
recently found in a more sophisticated calculation within 
the STLS approximation^! 

Conclusions - In summary, we have found the spin 
resolution of the electron gas correlation energy, via an 
approach applied to but not restricted to the three- 
dimensional uniform electron gas. Our results can be 
used to understand correlation in more realistic systems, 
and to construct improved density functionals and pair 
distribution functions. We have found that the same- 
spin correlation energy can be unexpectedly but under- 
standably positive. We have also provided support for 
resolutions^^ of the kinetic energy of correlation into f 
and J, terms. It is further possible to show that the pos- 
itive spin stiffness of correlation^ has positive tl and 
negative ft + II contributions. 
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